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Timing Signatures of the Internal- Shock Model for Blazars 

M. Bottcher^ and C. D. Dermer^ 

ABSTRACT 

We investigate the spectral and timing signatures of the internal-shock model 
for blazars. For this purpose, we develop a semi-analytical model for the time- 
dependent radiative output from internal shocks arising from colliding relativistic 
shells in a blazar jet. The emission through synchrotron and synchrotron-self 
Compton (SSC) radiation as well as Comptonization of an isotropic external ra- 
diation field are taken into account. We evaluate the discrete correlation function 
(DCF) of the model light curves in order to evaluate features of photon-energy 
dependent time lags and the quality of the correlation, represented by the peak 
value of the DCF. The almost completely analytic nature of our approach allows 
us to study in detail the influence of various model parameters on the resulting 
spectral and timing features. This paper focuses on a range of parameters in 
which the 7-ray production is dominated by Comptonization of external radia- 
tion, most likely appropriate for 7-ray bright flat-spectrum radio quasars (FS- 
RQs) or low- frequency peaked BL Lac objects (LBLs). In most cases relevant for 
FSRQs and LBLs, the variability of the optical emission is highly correlated with 
the X-ray and high-energy (HE: 100 MeV) 7-ray emission. Our baseline model 
predicts a lead of the optical variability with respect to the higher-energy bands 
by 1 - 2 hours and of the HE 7-rays before the X-rays by about 1 hour. We show 
that variations of certain parameters may lead to changing signs of inter-band 
time lags, potentially explaining the lack of persistent trends of time lags in most 
blazars. 

Subject headings: galaxies: active — gamma-rays: theory — radiation mecha- 
nisms: non-thermal 
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Introduction 



Blazars, a class of active galactic nuclei (AGNs) comprised of Flat- Spectrum Radio 
Quasars (FSRQs) and BL Lac objects, exhibit some of the most violent high-energy phenom- 
ena observed in AGNs to date. Their spectral energy distributions (SEDs) are characterized 
by non-thermal continuum spectra with a broad low-frequency component in the radio - 
UV or X-ray frequency range and a high-frequency component from X-rays to 7-rays. They 
show rapid variability across the electromagnetic spectrum. In extreme cases, the very-high- 
energy (VHE ) 7-ray emission o f blazars has been obse rved to vary on time scales of just a 
few minutes ( 1 Albert et al.l 120071 : lAharonian et al.l 120071 ) . 



The flux variability of blazars is often accompanied by spectral changes. Typically, the 
flux is most rapidly variable at the high-frequency ends of the two broad spectral components 
of the blazar SED. In the case of quasars, this refers to the optical (B-band) to UV and MeV 
to GeV 7-ray bands, while in the case of high-frequency peaked BL Lac objects (HBLs) it is 
the X-ray and VHE 7-ray regimes where the variability is the most extreme. In a few HBLs, 
the X-ray spectral variability could occasionally be characterized by s pectral hysteresis, i.e. 



a loop track of the blazar's state in a hardness-intensity dia gram (e.g. iTakahashi et al.lll996 



Kataoka et al.l l2000t iFossati et al.l I2OOOI : IZhang et al.l 120021 ) , although even within the same 
object this feature tends not to be persistent over multiple observations. Also in other types 
of blazars, hints of time lags betwe en different observing bands are occasi onally found in 
individual observing campaigns (e.g. 



Bottcher et al. 2007 



Horan et al.ll2009[ ) , but the search 



for time-l ag patterns persistin g throughout multiple years has so far remained unsuccessful 



(see, e.g. iHartman et al.ll200ll for a systematic search for time lags between optical. X-ray 
and 7-ray emission in the quasar 3C279). 

In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in 
a relativistic jet. The high-frequency (X-ray - 7-ray) emission could either be produced 
via Compton upscattering of low frequency radiation by the same electrons re sponsible 
for t he synchrotron emission (leptonic jet models; for a recent review see, e.g., iBottcher 
20071 ). or due to hadronic processes initiated by relativistic pr otons co-accelerated with 



the electrons (had ronic models, for a recent discussion see, e.g., iMiicke &: Protherod 12001 



Miicke et al.ll2003[ ). Leptonic models have been considered in a time- dependent manner with 



Kirk, Rieeer, & Mastichiadis 


1998; 


GeoreanoDoulos & Marschei 


1998; 


Chiaberee & Ghisellini 




1999; 


Kataoka et al. 


2000; 


Kusunose, Takahara & Li 


2000 


Li &; Kusunose 


2OO0I 


Bottcher & Chiang 


2002; 


Joshi & Bottcher 


2007 


. The time-dependent analysis of homogeneous single- zone lep- 



tonic models showed that spectral hysteresis patterns can be reproduced in a scenario of 
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gradual particle acceleration and subsequent radiative cooling, and that the presence and 
direction of hysteresis patterns depends on the relative values of the time scales for pa rticle 
acceleration, escape, and radiative cooling (jDermerlll998l : IChiaberge &: Ghisellinilll999[ ). 



Homogeneous leptonic jet models have met with remarkable success explaining the 
SEDs and correlated variability in many blazars. However, several recent observational 
results have seriously challenged homogeneous models and have motivated the consideration 
of inhomogeneous jet models. These observations i nclude the uncorrelated variability be- 
tween X-rays and 7-rays in the HBLs 1ES1959-I-650 (IKrawczynski et al.l 120041 ) and Mrk 421 
(IBlazejowski et al.ll2005[ ). and the uncorrelated opti cal and TeV emissi ons in PKS 2155-304, 
while X-rays and TeV 7-rays were well correlated ( ICostamantd l2008l ) . A particularly well 



motivated inhomogeneous blazar model is the internal shock model (e.g. ISpada et al.l 12001 



Sokolov. Marscher fc McHardvl l2004l : iMimica et al.l 120041 : iGraff et al.ll2008l ). In this model, 
the central engine is intermittently ejecting shells of relativistic plasma at varying speeds, 
which subsequently collide. Such models have had remarkable success in explaining SEDs 
and time lag features of generic blazars and deserve further study. 

The realistic treatment of radiation transfer in an internal-shock model for a blazar re- 
quires the time-dependent evaluation of retarded radiation fields originating from all parts of 
the shocked regions of the jet. The model system is th erefore highly non-linear and can gen- 



erally only be solve d using numerical s imulations ( e.g., ISokolov. Marscher fc McHardyll2004l : 
Mimica et al.l 120041 : iGraff et al.l 120081 : iJoshil 120091 ). As the current detailed internal-shock 



models employ either full expressions or accurate approximations to the full emissivities 
of synchrotron and Compton emission, a complete simulation of the time-dependent spec- 
tra and light curves is time-consuming and does therefore generally not allow to efficiently 
explore a large parameter space. General patterns of the SED, light curves and expected 
time lags between different wavelength bands have been demonstrated for very specific, but 
observationally very poorly constrained, sets of parameters. 

For this reason, we here develop a simplified internal-shock model in which the time- 
dependent synchrotron and external-Compton (EC) spectra are calculated completely an- 
alytically, and the SSC emission is reduced to a two-dimensional integral to be performed 
numerically. This approach allows us to calculate time-dependent snapshot spectra and light 
curves within a few minutes, and scan a large parameter space for the resulting spectral and 
timing features. 

Observational data of blazars often have limited, incomplete, and irregular time sam- 
pling, which complicates the evaluation of inter-band time lags and cross-correlations. The 
routinely used analysis m ethod designed to overco me these problems is the discrete corre- 
lation function (DCF, see lEdelson fc KrolikI Il988l ). Therefore, in order to produce results 
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directly comparable to observations, we subject our simulated light curves to the same DCF 
technique, and evaluate predicted inter-band time lags and the quality of the correlations, 
as represented by the peak values of the DCF. 

We describe the general outline of our model in ^ The dynamics of particle acceleration 
and cooling, and the resulting space- and time-dependent particle distributions, will be 
derived in ^ In §1] we describe our evaluation of the time-dependent radiative output from 
the internal-shock model. We present and discuss the results of a general parameter study 
in §3 and conclude in §H1 The appendices contain some details of the rather cumbersome 
analytical integrations required to evaluate the radiative output. 



2. Model Setup and Shell Dynamics 

We follow the collision of two relativistically moving shells (labeled a and b) in a blazar 
jet, powered by an intermittent central source. The basic geometry is illustrated in Fig. [T] 
(For a related trea tment of a relativistic shell interacting with a shell of material at rest. 



see 



Dermerl l2008l ) The two shells are being ejected from the central engine with Lorentz 
factors Ta,b with F;, > 1 and associated normalized velocities f3a,b (with (3 = v/c). 

In the rest-frame of the central engine, the ejection events of the two shells last for a time 
Ata fc. Consequently, the shells have widths (in the central-engine rest frame) of Ar^ f, = 
c(ia,b^ta,b- This assumcs that the shells do not spread along or transverse to the direction 



Demotion. The former effect can be important at ra,b ^ F^^cAta^f, (e.g., iMeszaros et al. 



19931 ). but is neglected because of the short duration of the colhsion during which the shell 
can be approximated as having a nearly constant thickness. Sideways expansion of the 
jet can also be neglected because of the short d uration of the co llision (this effect can be 



important for narrow decelerating jets; see e.g., ISari et al.l Il999l ). This latter assumption 
is also supported by observations of extragalactic jets remaining well colhmated out to kpc 
scales. If this assumption is vahd, the shell dynamics will not depend on the time between 
the shell ejections. 

Assuming that before the collision, the kinetic energy of the plasma in the shells' rest 
frames is negligible compared to the bulk kinetic energy ("cold plasma"), the kinetic lumi- 
nosity of the shells is given by = Fa feMa ^c^. This yields particle densities in the shells 

of ' ' ' L 

where R is the cross-sectional radius of the jet (and hence the shells), and the prime denotes 
quantities in the shell rest frame. 
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Fig. 1. — Schematic of the colhding-shell geometry. 

The colhsion between the two shells will lead to the formation of a forward shock moving 
into the slower shell a, and a reverse shock moving back into shell b. In the following, an 
overline denotes quantities in the frame of the shocked material behind the shock fronts. In 
this frame, the forward and reverse shocks move with Lorentz factors Tf r, respectively, away 
from the contact discontinuity. Pressure equilibrium across the contact discontinuity yields 
the condition 

n'a (r; - r^) = K (Tl - r.) . (2) 

Denoting the bulk Lorentz factor of the shocked material with respect to the stationary 
(AGN) rest frame as T (with velocity /3c), the shell velocities can be transformed to the 
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shocked fluid frame tlirougli 



ra,fer(i-/3,,,/3). 



(3) 



Eqs. 12] and 



Sari fc Piran 



can be solved to find T, and F 



1995: Panaitescu fc Meszaros 1999 



The solution is obtained numerically (cf. 



Dermdl2Q08|). 



The widths of the shells in the shocked-fluid frame are given by 

Ara,b = = 

J- f,r 



(4) 



where we have set /3afi = 1 since F^ b ^ 1. The time it takes for the shocks to cross their 
respective shells, can then be calculated as 



At 



a,b '-^l^a,b 



(5) 



in the shocked-fiuid frame. 



3. Relativistic Particle Dynamics 

Relativistic particles are entering the shocked-fiuid region with collective Lorentz factors 
Fj r, respectively. Upon shock crossing, the plasma will be compressed by a compression ratio 
r. For the results presented in ^ we have used r = 4. We assume that a fraction of 
the energy density in the shocked plasma will be contained in the magnetic field behind the 
shock fronts, which yields 



Bf,r = y Sttt eb (t)^^ - F/,,) < mpc2 (6) 

for the magnetic field strength behind the shocks. 

First- and second-order Fermi processes will accelerate particles at the shock fronts. 
We characterize the resulting injection of relativistic particles into the shocked-fiuid frame 
through the parameter eg giving the fraction of the shocked plasma kinetic energy in the 
shocked-fiuid frame that is transferred into relativistic electrons. Thus the injection power 
in relativistic electrons at the shock fronts is given by 



dt 



= ee irR^ mpC^ n'^^^ F/,^/3 (F/,^ - l) . (7) 

rel.e 



The injection of relativistic electrons is described as a power-law in electron energy, Ee 

■^f^ = Qil) = Qo 7"'^ ^(7; 71, 72) 
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where the triple-argument Heaviside function is defined as H{x] a,b) = 1 if a < x < b and 
otherwise, and all 7's refer to the shocked-fluid frame. We parameterize the width of the 
acceleration zone within the shocked-fiuid region as a multiple Aacc of the Larmor radius of 
a proton with 7p = Tf^r in the shocked-fiuid frame, 

Af acc;f,r = Aacc ^^^^^ ~ 3 X 10^ A^cc (^^^ T f,r CIR. (9) 

Consequently, the particle injection at any given point in the shocked-fiuid region will be 
active for a time 



Atacc;f,r = A,,, IllLH^f. ^ 10-4 A,,, (^) ' Tfl S. (10) 

Particle acceleration will hence be active in a volume V^acc;f,r = ^ri?^ Aracc;f,r- This allows us 
to evaluate the normalization of the particle injection function as 



- - - ( .4"'.-, if iq 7^ 2) 

VO;/,r - = X < (iij 

V/acc;f,r"^e -^L^- if g = 2. 



In pa 

VTl , 



The high-energy cutoff can be obtained by balancing the acceleration time scale of 
electrons with the synchrotron loss time scale. Writing the acceleration time scale as a 
factor Oacc times the electron gyration time scale, tacc = 27raacc7^eC^/(e-B), the maximum 
electron Lorentz factor will be given by 



72;/,, 



3e 



ttacc CTtB 



4.6 X 10^ a-}J^ 



B 



G 



-1/2 



:i2i 



If a fraction (e of electrons behind the shock fronts is accelerated into the power-law 
distribution ([8]) and we assume 72 ^ 7i, then the low-energy cutoff of the injection function 
([8]) is given by 



' mp q-2 e 
me q—1 Ce 



7l;/,r 



ff g > 2 



sst[r/.-i]72r/.r"^ ifi<g<2 



^|(r,,-l)/ln 



72;/,, 



if q 
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(13) 



where the q = 2 case will be solved numerically. 

The electrons injected according to Eq. [Slwill subsequently cool primarily due to radia- 
tive losses. Adiabatic losses become important after the shocks have traversed the shell and 
the shocked fluid shell begins to expand. Wh en this happen s, the magnetic field and there- 
fore the synchrotron emission rapidly decays (jDermerll2008l ). We neglect adiabatic losses in 
this study because of the rapid radiative cooling of the electrons. A good approximation to 
the time- dependent shape of the electron distribution can be found if all radiative losses can 
be described by a loss term 



7 = -^oT (14) 

which holds for synchrotron emission in a constant magnetic field as well as Compton scat- 
tering in the Thomson regime in a radiation field of constant energy density. In that case, 

i'Q = -caT ^ (15) 

where u is the sum of the magnetic-field and radiation field energy densities in the shocked- 
fluid frame. If the system under consideration is (a) SSC dominated and/or (b) Compton 
scattering at any given electron energy is predominantly happening in the Klein-Nishina 
regime, the evolution of the electron distribution can only be solved numerically. This 
is because in case (a), the dominant radiation field for Compton cooling depends on the 
current ( and recent) electron energy distribution and thus the cooling becomes non-linear 



see, e.g., iRoken fc Schlickeiserl l2009t ISchlickeiserl |2009| . for an analytical treatment of non 



linear radiative cooling), and in case (b), the cooling curve flattens towards higher energies 
compared to the simple 7^ dependence in Eq. [TH Therefore, our analysis is most directly 
applicable to the FSRQs and LBLs in which the 7-ray emission is generally believed to be 
dominated by EC scattering, mostly in the Thomson regime. 

If the conditions of a constant energy density u and Compton scattering in the Thomson 
regime are approximately fulfilled, then the time-dependent electron distribution can be 
found analytically. For this purpose, we define a spatial coordinate x which is defined in the 
shocked-fluid frame as x = at the contact discontinuity, and is measured positive in the 
direction of the forward-shock propagation (i.e., into shell a; see Fig. [1]). At any given point 
X, the time elapsed since the onset of the acceleration (i.e., since the forward or reverse 
shock has passed this point), is given by 
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tx^t-^ (16) 

where t denotes the reference time in the shocked-fluid frame. At any given point, the 
acceleration will remain active for a time Atacc;f,r- Then, the time- and space-dependent 
relativistic, non-thermal electron distribution can be well approximated by 

n{T, ^, tx) = Qo H{t^) -H'(7up - 7) |min(i^, At^^-i^^) 7"^ H{r, 71, 7c) 

+ Hi^ - 71) H{r, 7c, 7up) + At..e;f. 7r^ 7-^^(7; 7mi„, 7i) I (17) 

where the single- argument Heaviside function is defined as H{x) = 1 if a; > and otherwise, 
and the characteristic electron energies are given by 



1 



7up 



72 ^ + i^O max{0, {t^ - Aiacc;f,r)} 

1 



7c T 5 

l^O'f^x 

1 



7min _i J 

7l + ^otx 

(18) 



4. Radiative Output 

The evaluation of the radiative output from the entire shock structure will involve an 
integral of the emissivity along the x direction: 

uF^{e,tohs) = — ^2 — J ^h{x,t^,em) dx (19) 

^min 

where D = {T [1 — f3fj])~^ is the Dopplcr boosting factor with = cos6, the cosine of the 
viewing angle between the jet axis and the line of sight, and e = e{l + z)/D. The emission 
time tx,em has to be evaluated accounting for the light-travel time difference between different 
parts of the jet. The integration at a given observer's time tohs has to be performed in a way 
that dtohs — {1 + z) {dt — dx 11/ c) — 0, where t and x are measured in the stationary AGN 
frame. With di — dt/V and dx — Vdx, we find that an advancement in dx corresponds to 
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an advancement in emission time as dt = dx fi/{r'^c). Consequently, the co- moving emission 
time elapsed since the shock-crossing at any given point x at observed time corresponding 
to co-moving time t will be 

tx,em = t- + (20) 

The limits Xmin max of the integration in Eq. [12] will be given by constraints on the 
emission time being > (i.e., the shock front must have passed the point x), and the 
thickness of the shell: 



Xmin = - mm 



1 + 



2 



/ t(3fC > 

Xmax = niin = , Avn (21) 

\^l-/i/3^/r2 ' ^ ' 

There will be the additional constraint of particles being available to contribute to 
the emission at a given energy e, which depends on the individual emission mechanisms 
considered below. 



4.1. Synchrotron Emission 

For the purpose of our analytical treatment, we use a simple (5-function approximation 
for the synchrotron emissivity: 

where b = B/Bcnt with i?crit = {jrieC^) / [eh) ~ 4.4 x 10^^ G, and 7sy = a/c/^. Inserting this 
into Eq. [19] yields 



uF^yie, Us) = ^= ,J . e'/' / \ , U,em dx (23) 



v/TT^ 48 7r463/2 J My 5^ 

^min 

where we use Eq. [17] for the space- and time-dependent particle distribution ne(7sy, tx,em)- 
The integral in [23] can be solved fully analytically, and the solution is derived in Appendix 
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In this study, we focus on predictions for optical and higher-frequency emission. For this 
reason, we neglect synchrotron-self absorption (SSA) in our analysis. The emission region 
of our model system becomes optically thick to SSA at frequencies well below the optical 
regime. In our model systems considered in §3, the SSA frequency (where tssa = 1) is 
typically < ussa ~ 10^^ Hz. 



4.2. External-Compton Emission 

For the purpose of an analytical treatment, we evaluate the Comptonization of exter- 
nal radiation also with a simple (5-function approximation for the Thomson cross section. 
Furthermore, we assume that the external radiation field is isotropic in the stationary AGN 
frame and characterize it as mono-energetic with frequency i^ext, corresponding to a dimen- 
sionless photon energy in the co-moving frame, = F /iz/ext/('^eC^)- The radiation energy 
density Uext in the stationary AGN frame will be boosted into the shocked-fluid frame as 
Wext = r^Wext- The effect of the Klein- Nishina decline of the Compton cross section is ap- 
proximated as a hard cutoff in the scattered photon spectrum at ec = (-D/[l + zjls)- The 
beaming patterns of the intrinsically isotropic synchrotron and the external- Compton emis- 
sions a re slightly diff erent because the external radiation field is anisotropic in the co-moving 



frame (lDermerlll995l ). However, as long as the observer is located within the beaming cone 
at ^obs ~ 1/r, the difference is small and will be neglected in our simplified semi-analytical 
treatment. With these approximations, the treatment of the external Compton radiation 
is completely analogous to the one of synchrotron emission (see Appendix |A]), with the 
substitutions B"^/ (Stt) — )■ Wrad and b — )■ e^: 



^min \ / 

(24) 

with the solution to the integral given by the sum of the terms Ur + hf {i = 1,2, 3) for 
7 = ^y{e{l + z)/\es D]) derived in Appendix |Al 



4.3. Synchrotron Self-Compton Emission 



For the evaluation of the synchrotron-self-Compton emissivity je,sscy to use in Eq. [191 
we adopt, again, a 5-function approximation for the Compton cross section. 



- 12 - 



^ ar <5(e, - 7'[1 - /^/ije.) ^(^^c - ^^e) H{1 - 7e,[l - /J/lJ). (25) 

where /I^ is the cosine of the colhsion angle between the scattering electron and the incoming 
soft photon with energy e^. In Eq. |25]the effect of the Klein- Nishina decline of the Compton 
cross section is incorporated as a hard cutoff at the transition from the Thomson to the 
Klein-Nishina regime. With these approximations, the Compton emissivity becomes 



^ 1/(?[1-Mc]) 



J?,SSc(a;, tx,em) ~ — e ' \ ails / rfe^ V 1 - fl^ ^ ^e(7c, X, t^^em) 

(26) 



47r 

with 



For relativistic electrons, 7 ^ 1 we can neglect the factor /3 in Eq. [27] as long as 
/Ic ^ /3 ~ 1, so that it is an explicit solution. In order to use Eq. [261 we need an expression 
for the synchrotron photon density nph(e<i, f2s), which is a convolution of the (retarded) 
contributions from all shocked parts of the jet. For simplicity, we assume that all synchrotron 
photons enter a given point along the jet either directly from the front (superscript '+') or 
from the back (superscript '-'). In those cases, the scattering angle Ji^ is determined through 
the angular 6 distribution in Eq. [25] as 

where /i is the cosine of the observing angle in the observer's frame and /3r = a/1 — l/F^. 
We then write the synchrotron photon density as 



^ph(fs) ^x,cm) — '^r)ii(£s? ^x,em) „ ~l~ '^□^(f s; ^x,em) „ (^9) 

^ Zn ^ Zn 



where now the evaluation of ^ph(^s) involves an integral over the retarded emission from the 
shocked parts of the jet in front of and behind the point under consideration, respectively. 
This is requires the evaluation of an integral 
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■^s,max i? 



1 /" f T" diT 



s,min 



^s,max 

. (- -/^2 I d2 



^^'V(^')ln(^^7=^^=,^2^) (30) 



Sir c J \ {x — x') 



s,min 



The physically relevant case at hand here will correspond to thin slabs with R ^ Ara,b, for 
which the integral simplifies to 



•^s.max 



An c J \\x — X 



R 



^thi^) = I dx'hpY,{x') In ( 1^ ^^^,1 ) (31) 



The photon density distributions n^^ are evaluated fully analytically, as discussed in 
Appendix |Bl while the remaining two integrations over Is and x will be done numerically. 
Those are the only numerical integrations needed in our evaluation of the time-dependent 
emission spectra and light curves from the internal shock model. 

In our analysis, we neglect second-order SSC emission. This is justified by (a) the small 
Thomson depth (~ 10""^) of our model systems considered in the following section, and (b) 
the Klein- Nishina suppression of higher-order SSC emission. 



5. Results 

We have applied the semi-analytical internal-shock model described in the previous 
section, to parameter sets typical for FSRQs and LBLs. In these objects, the 7-ray emission 
is generally believed to be dominated by Compton scattering of external photons in the 
Thomson regime, in which case our approximate description of a constant electron cooling 
coefficient is valid. Table [1] lists parameters which we adopt for a baseline model in this 
regime. The chosen model parameters result in a bulk Lorentz factor of the shocked fiuid 
of r = 18.2, and a relative Lorentz factor between the shells of Frei = 1.133. The observing 
angle 6'obs has been chosen to coincide with the angle for which D = T, and we take the eg 
and es parameters equal in the forward and reverse shocked regions. 

Fig. |2] shows the instantaneous broadband spectra from this baseline model. In order 
to compare the SEDs with observations typically requiring extended exposure times, we also 
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evaluate the average SED over an integration time of 30 ksec. This time is representative 
of a deep blazar observation in X-rays, though shorter than typically required for Fermi to 
accumulate a signal with useful photon statistics from a bright 7-ray blazar. For Fermi, the 
minimum observation time from signal-dominated statistics for the detection of 5A^5 photons 
at 1 GeV photon energy is > 25A^5/(^i/5 ■ 8000 cm^ • 10"^° /-lo erg cm'^ s'^ ■ 1.6 x lO^^ 
ergs/GeV) ^ 5 x 10^A'"5/(Xi/5/_io) s, where O.2X1/5 is the fraction of time t hat Fermi will 



effect ively observe any given target in scanning mode covering about 2.4 sr (lAtwood et al. 



20091 ) . For the bright flare parameters used in Fig. 2, the uFi, flux at 1 GeV reaches values 
of ~ 5 X 10^'^ Jy Hz or /_io ~ 5, so that variability as short as ~ 3 hours could be detected. 
For comparison, the brightest blazar flare yet detected with Fermi was from 3C 4 54.3, which 



reached flux levels >2xl0 ^ph(>100 MeV) cm ^ s ^ (lEscande Sz Tanakall2009l ) . somewhat 
dimmer than the example considered here. 

The time-averaged SED is illustrated by the heavy solid curve in Fig. [2l It displays a 
synchrotron peak in the infrared, at ~ 2 x 10^^ Hz (A ~ 15/im), as well as an EC dominated 
7-ray peak at ~ 10^^ Hz (~ 40 MeV), with the 7-ray flux dominating over the synchrotron 
peak flux. It also displays a very hard, SSC-dominated X-ray spectrum [ax < 1)- These 
SEP properties are repr esentative of blazars of the FSR Q subclass, e.g., PKS 1510-089 



flD'Ammando et al.ll2009[ ) or 3C 454.3 flAbdo et al.ll2009bf ). We note that due to the hard 



low-energy cutoff introduced by our delta function approximation to the Compton emissivity, 
the emission in the hard X-ray - soft 7-ray (~ 100 keV - a few MeV) region may be 
underproduced in our simulations. 

The snap-shot SEDs in Fig. [2] illustrate that the high-energy end of the synchrotron 
peak, emitted directly behind the forward and reverse shock fronts, remains essentially un- 
affected as long as the observer receives synchrotron emission from the shocks still being 
located within the shells. As the shocks propagate, an increasingly larger region of the shells 
is energized with particles having longer time to cool. Therefore, the synchrotron spectrum 
extends to progressively lower frequencies. As the observer sees the shock regions leaving the 
shells, the highest-energy electrons rapidly cool off, leading to a loss of the high-frequency 
synchrotron emission. In the SSC component, the light-travel time delay leads to a delayed 
response of the SSC emission with respect to the synchrotron emission, with slightly cooled 
electrons still being able to efficiently Thomson scatter synchrotron photons up to 7-ray 
energies. 

Fig. [3] shows the light curves in the optical (R-band), X-rays (1 keV), and high-energy 
(HE) 7-rays in the Fermi range (100 MeV and 1 GeV) resulting from our baseline model sim- 
ulation. The synchrotron-dominated optical and the EC-dominated HE 7-ray light curves 
exhibit similar shapes, dominated by a rapid onset of the emission, a continued gradual 
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build-up as long as the shocks are located within the shells, and a rapid decay dominated by 
radiative cooling. The time scale of this decay is inversely proportional to the characteristic 
electron energy responsible for the respective emission. This explains the more rapid decay 
of the 1 GeV light curve compared to the 100 MeV one. The X-ray light curve, dominated 
by SSC emission, exhibits a much more gradual decay due to the broad-band convolution 
of electron and synchrotron photon energies involved to produce emission at any given fre- 
quency. At ~ 1.6 X 10^ s, the optical emission begins to be dominated by low-frequency SSC 
emission characterized by a slow, gradual decay. 

We point out that generally, blazars exhibit a low intensity of quiescent emission through- 
out the electromagnetic spectrum, which will prevent the observed light curves to reach the 
very low emission levels predicted for the beginning and very late times of our simulation. 
Assuming that this quiescent level of emission exhibits only moderate variability, it will lead 
to more moderate variability amplitudes, but will not affect the conclusions about cross 
correlations and time lags presented in the following. 

Even though some features of the SEDs and light curves in Figures [2] and E] are artifacts 
of our simplified analytical treatment of the time- and space-dependent particle distributions 
and the radiation processes, our results are expected to capture the salient spectral and light 
curve features as they would result from a more detailed numerical treatment. In particular, 
comparison is to be made with cross- correlation and time lag features that can realistically 
be extr acted from uneven l y sarn pled observational data. This is routinely done using the 



DCF of lEdelson fc KrolikI (|1988[ ). Thefore, for direct comparison with observational results, 
we apply the same formalism to our simulated light curves. The DCFs between the optical. 
X-ray and 100 MeV light curves shown in Figure [3l are plotted in Figure HI 

Figure S] illustrates that the more gradual decay of the (SSC-dominated) X-ray light 
curve with respect to the optical one translates into an optical lead before the X-rays by 
about 2 hr. The 100 MeV (EC) emission is produced by electrons of lower energy than 
the R-band optical synchrotron emission. This leads to a more gradual light curve decay 
at 100 MeV, translating into a 7-ray lag of ~ 1 hr behind the optical, and a reduced lag 
of X-rays behind the 100 MeV 7-rays (compared to the X-ray - optical lag) by only about 
~ 1 hr. All bands show a strong correlation with a DCF peak > 0.9. 

The characteristics of the injected electron distribution are evaluated according to Eqs. 
( IT3|) and (fT2|) . In our baseline model, these are 71 j = 130, ji^r = 350, and 72;/,r = 4.6 x 10"^. 
The magnetic field, according to Eq. Elis 1.0 G in both shocked regions. We note that the 
inferred va lues of 71 are substari tially higher than the c haracteristic values of 71 ~ a few 



inferred by iGhisellini et al.l (|2009[ ) and ISikora et al. J2009I ). In principle, similar values could 



also be achieved with different parameter choices in our simulations. However, in order to 
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still produce a synchrotron peak in the same range as observed in FSRQs, a substantially 
higher magnetic field {B > 10 G) would be required. Therefore, if the model explored here 
is applicable to FSRQs and LBLs with parameters close to the ones chosen in our parameter 
study, one would infer a substantially s maller total number of n onth ermal electrons and hence 



a smaller total jet power than found in iGhisellini et al.l ( l2009l ) and ISikora et al.l (120091 ) . 



Starting from our baseline model discussed above, we are now investigating the infiuence 
of the various parameters listed in Table [T] on the SED and light curve features. For each 
set of parameters, we evaluate the time-averaged SED over an integration time of 30 ksec to 
find the location of the synchrotron and Compton peaks as well as the Compton dominance, 
defined as CD = vFl^ / vF^^ . For each simulation, we calculate the DCFs between the optical 
(R-band), X-ray (1 keV), and HE 7-ray (100 MeV) light curves, and find the predicted time 
delay and peak value of the DCF. 

Specifically, we explore variations of the following parameters and their infiuence on the 
resulting time-averaged SED, inter-band time lags and DCF peaks: 

• The external radiation energy density, Wext 

• The electron equipartition fraction, eg 

• The magnetic-field equipartition fraction, eb 

• The fraction of electrons accelerated, Ce 

• The electron acceleration time scale parameter, aacc 

• The electron injection spectral index, q 

• The time scales of relativistic shell ejection, Ata,b 

• The relative Lorentz factors of the colhding shells, Frei ~ (l/2)(ra/r;, + Fb/Fa) 

• The kinetic power of the faster shell, Lf, 

• The radius of the shells, R 



Each parameter is varied individually, leaving all other parameters constant at the value 
of our baseline model (Tabled]). In the following, we discuss the infiuence of each individual 
parameter mentioned above, on the SED and DCF characteristics. 
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5.1. Variations of Uext 



An increasing external radiation energy density leads to a transition from SSC to EC 
dominated 7-ray production and overall more rapid electron cooling. The consequences 
in the SED are decreasing synchrotron and EC peak frequencies, along with a decreasing 
synchrotron pe ak flux and an incre asing Compton peak flux (see Fig. E]). This shift in SED 
characteristics (IFossati et al.lll998l ) reflects the gradual tra n sition from BL-Lacs to FSRQs 
along the blazar sequence as suggested by iGhisellini et al.l ( 1l998l ). The dependence of the 
external fleld energy density shown in Fig. [5]would only apply to the FSRQs and LBLs, which 
have peak synchrotron frequencies between 10^^ and 10^^ Hz. In order to extend this to XBLs 
and the blazar sequence, it may also be n ecessary to posit a correlation between jet power 
and the external broad-line region density (iBottcher fc Dermerll2002l ; IGhisellini &: Tavecchio 
2008h . 



Fig. [H] illustrates the dependence of time lags and DCF peak values on the external 
radiation energy density. For high values of Ue^t 3 x 10""^ erg cm~^), the more rapid 
cooling with incresing Uext leads to generally decreasing time lags. For smaller values of 
Mcxt, we begin to see the effects of SSC emission becoming signiflcant in the 7-ray regime. 
Most notably, this ultimately leads to an inversion of the optical - HE 7-ray lead into a lag 
of optical behind HE 7-rays for very low Wext- The reason for this is that the SSC 7-rays 
at 100 MeV are produced by electrons of higher energy than those emitting synchrotron 
radiation in the R band. 

All three bands remain generally well correlated (DCF peak > 0.75), with only a weak 
trend towards lower correlation quality for SSC dominated high-energy emission. 



5.2. Variations of eg 

With increasing electron equipartition fraction eg, the total energy input into relativistic 
electrons increases, along with an increase in the low-energy cutoff of the injected electron 
distribution. The increased electron energy density leads to an increasing radiative output 
in both spectral components (synchrotron and Compton) as well as an increasing fraction of 
SSC to synchrotron emission. Since the HE 7-ray emission remains EC dominated for our 
parameters, the Compton dominance remains essentially unchanged. The peak frequencies 
of the two spectral components remain essentially unchanged for realistic values of eg ^ 0.1 
(see Figure [7]). 

There is an overall weak trend of decreasing inter-band time lags with increasing eg, 
which is related to decreasing electron cooling time scales. All three bands remain well 
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correlated (DCF peak > 0.85) irrespective of eg. 

5.3. Variations of 

The parameter regulates the magnetic field. An increasing value of implies stronger 
synchrotron cooling. In the synchrotron frequency range, this effect is largely cancelled out 
by the increase in the characteristic synchrotron frequency {usy oc Therefore, the 

synchrotron peak frequency remains essentially unchanged as changes, while the Compton 
peak frequency decreases (see Fig. |9]). 

As a consequence of the limited amount of power available to channel into synchrotron 
and Compton emission, the Compton peak flux remains almost constant with varying es, 
while the synchrotron peak flux continuously increases, consequently leading to a decreasing 
Compton dominance. 

In the DCFs (see Figure [TOj) . there is a general trend of decreasing time lags and in- 
creasing quality of correlations with increasing e^. This is a consequence of the decreasing 
radiative cooling time scales (with progressively larger contributions from synchrotron and 
SSC cooling). 

5.4. Variations of 

Increasing the fraction (e of electrons accelerated, while keeping the power injected into 
relativistic electron (parametrized by eg) constant, leads to a decreasing low-energy cutoff 
of the electron distribution. As Ce increases, an increasing fraction of the available power 
in electrons will be distributed at low electron energies, leading to a decreasing radiative 
output. Consequently, both the synchrotron and Compton peak frequencies as well as the 
synchrotron and Compton peak fluxes decrease with increasing (e. 

Less efficient radiative cooling because of lower average particle energies leads to a 
gradual trend of increasing time lags with increasing (g. 

5.5. Variations of Oacc 

The electron acceleration efficiency parameter Oacc regulates the high-energy cutoff of 
the injected electron population via Eq. [121 A lower value of Oacc implies a larger value 
of 72. Changes in Oacc (and hence 72) have a negligible influence on the location of the 
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sychrotron and Compton peaks, both in frequency and flux (for electron spectral indices 
q > 2). However, the value of Oacc (and hence, 72) determines whether the synchrotron 
emission of freshly injected electrons extends into the X-ray regime or cuts off at optical- 
UV frequencies. As Oacc increases, the synchrotron emission cuts off at progressively lower 
frequencies. Specifically, this leads to an increasingly rapid decline of the R-band emission 
after the shocks have broken out of the shells. As long as the synchrotron emission of freshly 
emitted electrons does extend beyond the R-band, this leads to an increasing R-band lead 
before the X-ray emission. For a very high value of Oacc (10^ in our simulations), even the 
synchrotron emission of freshly injected electrons cuts off at frequencies close the R-band, 
and the R-band flux is dominated by low-frequency SSC emission very early on. This reverses 
the R-band vs. X-ray lead into a lag. A similar argument applies to the dependence of the 
100 MeV 7-ray emission with respect to X-rays and optical. 

All light curves correlate well with DCF peakds > 0.9, with no substantial dependence 

on Oacc- 

5.6. Variations of q 

The electron spectral index at the time of injection, q, obviously directly determines the 
spectral shape of the synchrotron and Compton emissions. In particular, for hard injection 
spectra {q ^2), the peak frequencies are determined primarily by the energy of the highest- 
energy electrons. As we choose a steeper injection electron spectrum, the peak shifts to a 
position where it is predominantly determined by the low-energy cut-off 71 and thus only 
weakly dependent on the electron injection index (see Figure [T^ . 

Generally, a harder electron injection spectrum (lower q) implies that a larger fraction 
of the energy transferred to electrons is stored in high-energy electrons. This leads to a 
monotonic increase of the Compton dominance with decreasing q (harder electron spectrum). 

Generally, all frequency bands correlate well with DCF peak values > 0.8. There is a 
slight trend of de-correlation of the HE 7-rays with optical and X-rays for hard injection 
spectra {q < 2.3). The increasing quality of correlation with softening injection spectra goes 
in tandem with decreasing absolute values of the time lags. Notably, for very steep injection 
spectra (g > 3), the instantaneous synchrotron emission in the R band has a very steep 
spectral index so that a substantial fraction of the R-band flux is contributed by very-low- 
frequency SSC emission. This reverses the R-band lead observed for harder injection spectra 
into an R-band lag behind X-rays and HE 7-rays. 
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5.7. Variations of Ata.b 

The shell ejection time scale determines the width of the shells. This in turn, regulates 
the shock propagation time and thereby the dynamical time scale of the shocked emission 
region. Furthermore, since we left the ejection power constant, a longer ejection time corre- 
sponds to a larger total kinetic energy deposited into the shells. This latter effect causes the 
overall radiative power output to increase with Ata.b- The Compton peak is dominated by 
EC, so that the Compton dominance remains essentially unchanged. The values of the peak 
frequencies show only a weak increasing trend with increasing Ata,b- 

The most notable trend in the DCF features is a reversal of the lag of R behind HE 
7-rays (100 MeV) into a lead for long shock propagation time scales, in tandem with an 
increasing X-ray lag behind HE 7-rays. This is due to the fact that both X-rays and R-band 
are dominated, at least at late times, by slowly decaying SSC emission. As the shocked 
regions are more extended, SSC emission persists for a longer time due to light-travel time 
delays pertinent to the SSC emission. In contrast, the decay of the EC-dominated HE 7-ray 
emission is determined only by the radiative cooling time scale, which is independent of 
Ata,b with the assumptions made here. It is notable that the time lags are generally not 
proportional to the dynamical time of the shock propagation, as one could naively expect. 

5.8. Variations of Frei 

We have explored the effect of an increasing difference in the Lorentz factors of the 
colhding shells, leading to an increasing relative Lorentz factor, Frei ~ (1/2) (Ta/Tb + Tb/Ta), 
by varying Fq and F;, in a way that the resulting shocked-fluid Lorentz factor F (and hence 
the Doppler boosting factor D) remains unchanged. An increase in Frei drastically increases 
the radiative efficiency of the shells. The increased energy density in the shocked shells will, 
at the same time, increase the magnetic fields and decrease the radiative cooling time scales. 
It will also lead to larger shock propagation Lorentz factors Fj ,. and hence larger low-energy 
cut-offs in the electron distributions. This leads to a net increase of the peak frequencies of 
both radiation components with increasing Frei. The peak fluxes in both components increase 
drastically with increasing Fj-ei with only a moderate effect on the Compton dominance. 
The decreasing radiative cooling time scales with increasing Fj-ei imposes a general trend of 
shortening inter-band time lags with only minor effect on the DCF peak values (see Fig. [T5|l . 
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5.9. Variations of Lf, 

Varying the kinetic luminosity of one of the shells, keeping its Lorentz factor unchanged, 
directly amounts to increasing its particle density. For unchanged e^, this will also increase 
the magnetic field. Both effects lead to an increasing synchrotron peak flux and a decreasing 
Compton dominance. Due to more rapid radiative (synchrotron) cooling, both peak peak 
frequencies decrease with increasing Lj,. As the system becomes increasingly synchrotron + 
SSC dominated, the hght-travel time delays inherent to SSC emission lead to a general trend 
of increasing time delays with increasing Lf,, accompanied by very moderate changes in the 
DCF peak values. 

5.10. Variations of R 

Increasing the shell radius implies a smaller particle density and hence a lower magnetic 
field. Consequently, the synchrotron peak flux and peak frequency decrease with increasing 
shell radius, leading to an increasing Compton dominance. 

In the DCFs, an increasing radius it! leads to a gradually increasing trend of time lags 
between R and higher-frequency bands with shghtly decreasing DCF peak values (~ 0.95 — > 
0.8 for R — 3 X 10^^ cm — > 3 x 10^^ cm). The increasing time lag trend may be related to 
the increasing synchrotron cooling time as R increases. 

6. Summary and Conclusions 

We have developed a semi-analytical internal-shock model for blazars. Our model in- 
cludes synchrotron, SSC, and external Compton emission, and properly takes into account 
light-travel time and shock-propagation effects as well as the space- and time-dependent 
relativistic electron distributions in the shocked shell region to evaluate the time-dependent 
radiative output. For direct comparison with sparsely and unevenly sampled observational 
data, we evaluated the Discrete Correlation Functions between several representative wave- 
lengths, namely between the optical R band, 1 keV X-rays, and high-energy 7-rays (100 MeV 
— Fermi). 

Our calculations were used to characterize inter-band time lags and the quality of the 
correlation, as represented by the peak value of the DCF. For this calcuation, we considered 
a specific model blazar where spectral variability could be detected with optical and pointing 
X-ray telescopes at the derived flux level. These calculations can be extended to apply to 
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brighter blazars where variabihty might be detected with the Fermi telescope, or otherwise 
used to model the time-averaged SED observed with Fermi. Keeping these caveats in mind, 
we have studied the influence of variations of a number of essential model parameters on the 
resulting SED and DCF characteristics. 

We note that our model results apply to FSRQs and LBLs, maybe IBLs, but not HBLs, 
because of the restriction to the Thomson regime. The parameter choices in our model are 
extensive, including eg and parameters, shell widths, collision radii, and jet powers. We 
have reduced the allowed parameter space by taking equal shell powers, with the forward 
and reverse shocks having equal eg and es parameters. By examining the dependence of 
peak frequency, peak vF^ flux and Compt on dominance in F ig. [5l we flnd a trend in accord 
with the behavior of the blazar sequence (iFossati et al.lll998[ ). 



Our results do not directly explain the break in the vFy sp ectra at a few GeV in blazars 
like 3C 454.3 and AO 0235+164 and other FSRQs and LBLs kbdo et all l2009al jbl) . or the 
relative constancy of the GeV spectral index with flux state. The GeV breaks in several 
blazar SEDs can not be explained as a cooling break since, whenever photon statistics allow 
a reliable determination of the break, it is Aa > 0.5. However, a smooth high-energy 
cutoff could effectively mimic a spectral break, and such a cutoff around a few GeV could 
be produced in our model by a lower high-energy cutoff of the electron distribution, which 
could be achieved with a choice of a larger value for our paramter Oacc- 

We have not considered a detailed spectral comparison with the impressive multiwave- 
length SEDs now availab le due to campaigns involving the Fermi Gamma Ray Space Tele- 
scope ( lAbdo et al.ll2009d ). Rather, here we considered a parameter study of short timescale 
variability between optical. X-ray, and gamma-ray eneriges. 

In most cases, the optical light curves showed substantial leads (by ~ a few hours) 
before X-ray and HE 7-ray emission. However, variations of several parameters (e.g., the 
external radiation energy density Wext, the acceleration time scale parameter Oacc, the electron 
spectral index g, and the shell widths) can drastically change the amount and even sign of 
the time delays. This may explain the lack of a consistent pattern of inter-band time la gs 
when comparing multiple observing epochs of the same object (e.g.. 



Hartman et al.ll2001f ) 



MB acknowledges partial support from NASA through INTEGRAL Guest Investigator 
grant NNX09AI71G and Fermi Guest Investigator grant NNX09AT82G. The work of ODD 
is supported by the Office of Naval Research. 



- 23 - 



REFERENCES 

Abdo, A. A., et al. 2009, ApJ, 700, 597 
Abdo, A. A., et al. 2009, ApJ, 699, 817 
Abdo, A. A., et al. 2009, ApJ, 697, 934 
Aharonian, F., et al., 2007b, ApJ, 664, L71 
Albert, J., et al., 2007, ApJ, 669, 862 
Atwood, W. B., et al. 2009, ApJ, 697, 1071 
Blazejowski, M., et al., 2005, ApJ, 630, 130 
Bottcher, M., & Chiang, J., 2002, ApJ, 581, 127 
Bdttcher, M., & Dermer, C. D. 2002, ApJ, 564, 86 

Bottcher, M., 2007a, in proc. "The Multimessenger Approach to Gamma-Ray Sources", 
ApSS, 309, 95 

B6ttcher, M., et al, 2007, ApJ, 670, 968 

Chiaberge, M., & Ghisellini, C, 1999, MNRAS, 306, 551 

Costamante, L., 2008, IJMP D, 17, 1449 

D'Ammando, F., et al, 2009, A&A, in press f' arXiv:0909.3484ll 

Dermer, C. D., 1995, ApJ, 446, L63 

Dermer, C. D. 1998, ApJ, 501, L157 

Dermer, C. D. 2008, ApJ, 684, 430 

Edelson, R. A., & Krolik, J. H., 1988, ApJ, 333, 646 

Escande, L., & Tanaka, Y. T., 2009, Astronomer's Telegram 7^ 2328, 1 

Fossati, G., et al., 2000, ApJ, 541, 166 

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 1998, MNRAS, 299, 
433 

Georganopoulos, M., & Marscher, A. P., 1998, ApJ, 506, 621 



Ghisellini, G., Celotti, A., Fossati, G., Maraschi, L., & Comastri, A., 1998, MNRAS, 301, 
451 

Ghisellini, G., Tavecchio, F., & Ghirlanda, G. 2009, MNRAS, 399, 2041 
Ghisellini, G., & Tavecchio, F. 2008, MNRAS, 387, 1669 

Graff, R B., Georganopoulos, M., Perlman, E. S., & Kazanas, D., 2008, 689, 68 

Hartman, R. C., et al., 2001, ApJ, 558, 583 

Koran, D., et al., 2009, ApJ, 695, 596 

Joshi, M., & Bottcher, M., 2007, ApJ, 662, 884 

Joshi, M., 2009, Ph.D. Thesis, Ohio University 

Kataoka, J., et al., 2000, ApJ, 528, 243 

Kirk, J. G., Rieger, F. M., & Mastichiadis, A., 1998, A&A, 333, 452 

Krawczynski, H., et al., 2004, ApJ, 601, 151 

Kusunose, M., Takahara, F., & Li, H., 2000, ApJ, 536, 299 

Li, H., & Kusunose, M., 2000, ApJ, 536, 729 

MAGIC Collaboration, et al. 2008, Science, 320, 1752 

Meszaros, P., Laguna, P., & Rees, M. J. 1993, ApJ, 415, 181 

Mimica, P., Aloy, M. A., Miiller, E., & Brinkmann, W., 2004, A&A, 418, 947 

Miicke, A., & Protheroe, R. J., 2001, Astropart. Phys., 15, 121 

Miicke, A., et al., 2003, Astropart. Phys., 18, 593 

Panaitescu, A., & Meszaros, P. 1999, ApJ, 526, 707 

Roken, C, & Schlickeiser, R., 2009, A&A, 503, 309 

Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 519, L17 

Sari, R., & Piran, T. 1995, ApJ, 455, L143 

Schhckeiser, R., 2009, MNRAS, 398, 1483 



- 25 - 



Sikora, M., Stawarz, L., Moderski, R., Nalewajko, K., & Madejski, G. M. 2009, ApJ, 704, 38 

Sokolov, A., Marscher, A. P., & McHardy, I. A., 2004, ApJ, 613, 725 

Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A., 2001, MNRAS, 325, 1559 

Takahashi, T., et al., 1996, ApJ, 470, L89 

Zhang, Y. H., et al., 2002, ApJ, 572, 762 



A. Time-dependent Synchrotron Spectra 

Since the particle energy distribution of Eq. [T71 explicitly depends on tem, it is convenient 
to transform the integration fl2^ into an integration over emission time, using 

As mentioned above, the limits Xmin,max of Eq. ED will then be augmented by the condition 
that n(7sy,tcm) > for a given electron energy 7sy. Since /x/F^ < l//3/,r in the parameter 
range relevant to internal shocks in blazars (/i ~ 1; F >> 1, (3j^^ ~ 1), the emission times 
will be < t = tobs(l + throughout the integration region. Thus, defining a time limit 
tf corresponding to the limit x^ax in the forward shock region and corresponding to x^nm 
in the reverse shock region, the integral in Eq. [23] will consist of two branches, 




The limits on x in Eq. [2T] can then be converted into limits on tx,em as 



max 



mm 



tl3j.c 



1 + 



r2 



fx 1 



This preprint was prepared with the A AS lAT^TjX macros v5.2. 
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tf = max < , I mill 

{A3) 



tl3fC 

_ !^ 

r2 



/i 



We now evaluate each term in the integral (lA2p separately for each term in the expression 
for n in Eq. (flTll : 



= ^ / -f^(7up - 7) H{^; 7i,.r, 7c) min(t 

x,em5 '^^acc,r) '^^x,em* 

(A4) 



_i_ J. 



This integral contributes only as long as 



^x,em ^ tr {'^) 

t^,em < t (b) 

7 < 7c = — (c) 

7 < 7up — max(0,[tx,om-Alacc.r]) 

(A5) 

These constraints translate into effective limits of the integral f lA4p of 



and hence, 
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The forward-shock contribution of the first term is 
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t 



hf = 7 ^ / H{'y^p - 7) if (7; 7i,r, 7c) min(tx,em, Atacc.r) C?tx,em- (A8) 

Analogous to the previous integral, the effective limits of this integration are 



and hence, 
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^1/ = ^ H{l - 7i,f) 



j_ 



,max ^lf,min) 
(l/2)(AtL,f-^f,„.in) + A^acc,f (t 
AtacCjf ('^lf,max '^lf,min) 



if ^lf,min ^ ^lf,max ^ Atacc,f 
If, max Atacc,f) if ^lf,min ^ Atacc,f ^ ^lf,max (AlO) 



if AtacCjf ^ '^lf,min ^ '^lf,max 

The reverse-shock contribution from the second term in Eq. [T7] is 



J _QoC-i ^"^^^ f min(tx,em,Aiacc,r) ,T 

hr = — / ii (7; 7c, 7up) H{'j - 7i,,) = dtx,. 



(All) 



with the limits determined by the conditions 



^x,em ^ tr {'^) 

^x,em < t (b) 

7 > 7c = — (c) 

7 < 7up - ^-,'+!/0 max(0,[Ix,om-Atacc.r]) 

(A12) 

These constraints translate into effective limits of the integral (lAlip of 
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i2r,max = mill <! t 
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(A13) 



and hence, 



'2r = ^ ^ ^(7 - 7l,rJ 
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Analogously, for the forward-shock contribution from the second term in Eq. [T71 we 
have the effective limits 
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The third term in Eq. [T7] yields a reverse-shock contribution of 
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with the hmits determined by the conditions 
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These constraints translate into effective hmits of the integral flA17p of 
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and hence. 
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Finally, the contribution of the third term in Eq. [T7] from the forward shock is 

-'3f = 7 X l%,max ~ %,minj K^"^^) 

with the limits 
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B. Evaluation of the synchrotron photon density for SSC 

For each incoming photon direction (±), the seed space- and time-dependent seed pho- 
ton density n^^^{esX,tx,em) consists of contributions from all three branches of the electron 
distribution f lT7|) . and may contain contributions from both the forward and reverse shock. 
For any given branch of the electron distribution, the calculation is analogous for the various 
cases to be considered. Therefore, we give here only one representative case for each branch. 

The space- and time-dependent photon density distribution is evaluated as 



s.mm 



(Bl) 

where 



X' ^ / / 1 'y* I 'y*^ /y* I 

Although B{x') and b{x') are constant throughout the forward shock region and the reverse 
shock region, we need to take into account that generally Bf ^ B^. The integration limits 
in Eq. IB II are set by the conditions 
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and hence 
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Because of the absolute values involved in (Eq. IB2p . we need to evaluate the integral 
in Eq. IB II separately for the cases x > (A) and x < (B). 

For case (A), the integral in Eq. IB II can be split up into 6 contributions: 

If^ ^ + If,^ + If,^) (B6) 

corresponding to the forward and backward traveling photons from the three terms in the 
expression for the electron density in Eq. [T71 Hence, 



2 "^s.max 

^si^ = Q^J-^^'' j ^(7up-7)^(7;7i,7c) min{t,,Ata,e,f} ln(^^^) (ix' (B7) 



with 7 = y^ts/hf. 



In the region corresponding to the A+ contribution, we parameterize the synchrotron 
emission time as 



= «A+ - /5a+ x' (B8) 
with 

«A+= t+f_(^ + l) 

/5a. = ^ (B9) 

The effective limits of the integration ( IB 71) are determined by the Heaviside functions 

as 
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To evaluate the integral I^^ with these limits, we define a critical x'^ for which tx{x'^ 

A^acc,f • 
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with tx < Atacc.f if x' > Furthermore, we write 

Then, the integration yields: 
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Analogous calculations (carefully accounting for the different implications of the absolute 
values in Eq. IB8I and the different magnetic-field values in the forward- and reverse-shock 
regions) yields the contributions A- and B± from the first term in Eq. [T71 



The A+ {x' > X > 0) contribution of the second term in Eq. [TTlis 
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The effective limits of this integration are 
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Using the parametrization for tx from Eq. IB8I and the critical value of x' from Eq. IBllI for 
which tx = Ataccf, the solution to is 
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is the Lerch function. 
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Again, analogous calculations yield the A- and B± contributions of the second term in 
Eq. HZl 

The A+ {x' > X > 0) contribution of the third term in Eq. [17] is 
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The effective integration limits are given by 
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and the result of the integration is then 
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Analogous solutions are obtained for the A- and B± contributions. 
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Fig. 2. — Snap-shot SEDs for our baseline model with parameters listed in Table [H The 
heavy solid (red) curve shows the SED resulting from averaging over an integration time of 
30 ksec, representative of a typical exposure time of a deep X-ray observation of a blazar. 
Individual radiation components of the time-averaged SED are shown as: dotted = syn- 
chrotron; long-dashed = SSC; dot-dashed = EC. The dotted vertical lines indicate the fre- 
quencies (R-band, 1 keV, 1 MeV, 100 MeV) at which light curves have been extracted. 
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Table 1. Parameters of our baseline model. 



Pcir 3,111 Gtcr 


Symbol 




Lorentz factor of shell a 




15 


Lorentz factor of shell b 




25 


Kinetic power of shell a 


La 


10'*^ erg s"'^ 


Kinetic power of shell b 


Lb 


10''^ erg s"'^ 


Duration of ejection of shell a 


Ata 


2 X 10^ s 


Duration of ejection of shell b 


Aii, 


2 X 10^ s 


Observing angle 




3.15° 


Shell radius 


R 


3 X 10^^ cm 


Electron equipartition fraction 




0.1 


B-field equipartition fraction 




10~^ 


Electron acceleration fraction 


Ce 


0.01 


Acceleration time scale par. 


'2'acc 


10^ 


Acceleration length parameter 


^acc 


10 


Elect, injection spectral index 


1? 


2.6 




O.Oe+00 5.0e+03 l.Oe+04 1.5e+04 2.0e+04 2.5e+04 3.0e+04 3.5e+04 

Us 

Fig. 3. — Light curves at various frequencies/energies resulting from our baseline model with 
parameters listed in Table [TJ 
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Fig. 4. — Discrete correlation functions between all energy bands for which light curves are 
plotted in Figure O A positive delay of "band 1 vs. band 2" indicates a lead of band 1 
before band 2. 
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Fig. 5. — Characteristics of the time- averaged SED as a function of the external radiation 
energy density Wext- 
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Fig. 6. — Characteristics of the DCFs as a function of the external radiation energy density 
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Fig. 7. — Characteristics of the SEDs as a function of the electron equipartition fraction eg. 



-41 - 




I I I I I I I I I I I I 

0.01 0.1 



Fig. 8. — Characteristics of the DCFs as a function of the electron equipartition fraction e^- 
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Fig. 9. — Characteristics of the SEDs as a function of the magnetic- field equipartition fraction 
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Fig. 10. — Characteristics of the DCFs as a function of the magnetic-field equipartition 
fraction €3- 
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Fig. 11. — Characteristics of the DCFs as a function of the electron acceleration efficiency 
parameter aacc- 
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Characteristics of the SEDs as a function of the electron injection spectral index 
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Fig. 13. — Characteristics of the DCFs as a function of the electron injection spectral index 
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Fig. 14. — Characteristics of the DCFs as a function of the shell ejection time scales, Ata.b 
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Fig. 15. — Characteristics of the DCFs as a function of the relative Lorentz factor between 
the shells, Frei- 



